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Abstract 



The aim of the present work is to investigate tokamak equilibria with reversed magnetic 
shear and sheared flow, which may play a role in the formation of internal transport 
barriers (ITBs), within the framework of two-fluid model. The study is based on ex- 
act self-consistent solutions in cylindrical geometry by means of which the impact of the 
magnetic shear, s, and the "toroidal" (axial) and "poloidal" (azimuthal) ion velocity com- 
ponents, Vi Z and Vig, on the radial electric field, E r , its shear, \dE r /dr\, and the shear of 
the E x B velocity, c^exb = \d/dr(E x B/P 2 )|, is examined. For a wide parametric 
regime of experimental concern it turns out that the contributions of the V{ Z , Vie and 
pressure gradient (VP,) terms to E r , \E' r \ and u ExB are of the same order of magnitude. 
The contribution of the VPj term is missing in the framework of magnetohydrodynamics 
(MHD) [G. Poulipoulis et al. Plasma Phys. Control. Fusion 46 (2004) 639]. The impact 
of s on cl>exB through the VPj term is stronger than that through the velocity terms; in 
particular for B z = constant, the contributions of the VP and velocity terms to cuexb 
at the point where dE r /dr = are proportional to (1 — s)(2 — s) and (1 — s), respec- 
tively. The results indicate that, alike MHD, the magnetic shear and the sheared toroidal 
and poloidal velocities act synergetically in producing electric fields and therefore o>exb 
profiles compatible with ones observed in discharges with ITBs; owing to the VPj term, 
however, the impact of s on E r , \E' r \ and ujexb is stronger than that in MHD. 
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Introduction 



Tokamak discharges with improved energy and particle confinement properties in con- 
nection with internal transport barriers (ITBs) have certain attractive features, such as a 
large bootstrap current fraction, which suggest a potential route to steady-state mode of 
operation desirable for fusion power plants. Long quasi-steady or steady ITB states have 
been obtained in different tokamaks, e.g ASDEX Upgrate [UI2], JT-60U [3], Tore-Supra 
[3], and JET [3] where ITBs were maintained for up to 11 s. The ITBs usually are associ- 
ated with reversed magnetic shear profiles jE|,[Z| and their main characteristics are steep 
pressure profiles in the barrier region [E] and radial electric fields associated with sheared 
flows [01 ED]. The mechanism responsible for the formation of ITBs and the underlying 
physics is not completely understood. Most theoretical models supported by experimental 
observations rely on suppression of microinstability induced transport in connection with 
reversed magnetic shear, s < 0, sheared flow, the radial electric field, E r , its shear, \E' r \, 
and most importantly the E x B velocity shear, 



WExB 



d E x B 



dr B 2 



In particular the E x B velocity shear may lead to a reduction in the amplitude of 
turbulent fluctuations, even to their suppression, or to a decrease in the radial correlation 
lengths [IT] . Although there are experimental observations supporting this scenario, the 
overall experimental evidence up to date is rather complicated, not universal in the various 
tokamak machines and has not made clear whether the magnetic shear or the sheared flow 
(toroidal or poloidal) are more important for the ITB formation. A discussion on this 
issue is made in the Introduction of Ref. [T2j. Also, the experimental and theoretical 
knowledge on discharges with ITBs was reviewed recently in Refs. an d [T3] . 

In a previous work ^2] we studied magnetohydrodynamic (MHD) equilibrium states 
with reversed magnetic shear and sheared flow in cylindrical geometry. In particular, 
presuming that E r , E' r and o>exb are of relevance to the formation of ITBs we examined 
how these quantities are affected by the magnetic shear and sheared flow and found that 
the latter quantities act synergetically in increasing c^exb with the impact of the flow, in 
particular the poloidal one, being stronger than that of the magnetic shear, s. The present 
work aims at extending the study to the framework of the two fluid model. This model is 
advantageous over MHD in that the contribution of the ion pressure gradient (V-Pj) term 
to E r , contribution which is missing in MHD, can be obtained from the ion (or electron) 
momentum equation. In addition the current density can be expressed self-consistently in 
terms of the ion and electron fluid velocities. Also, we shall examine the impact of certain 
local characteristics of the safety factor profile, i.e. the minimum of q and its position, 
on the aforementioned quantities (not addressed in Ref. ^2]) and the relative sign of the 
"toroidal" (axial) ion velocity, v iz , "poloidal" (azimuthal) ion velocity, v^, and toroidal 
magnetic field, B z . It turns out that, owing to the VP« term, the impact of s oia E r , E' r and 
WexB is stronger than that in MHD. The contribution of both flow components, however, 
remains significant. In addition in many cases s enhances the velocity contribution to 
these quantities which, alike in MHD, indicates a synergism of s and the flow. 

The work will be conducted through the following steps. Exact solutions of a slightly 
reduced set of two-fluid equilibrium equations for a cylindrical magnetically confined 
plasma are constructed in section 2 by prescribing the profiles of certain free quanti- 
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ties, including the safety factor and the toroidal and poloidal ion velocities, in accord 
with ITB experimental ones. Then in section 3 we examine the impact of s, the velocity, 
the velocity shear, the local characteristics of q, and the relative signs of the velocity 
components and B z on E r , E' r and o;e x b- The characteristics of the pressure and toroidal 
current density are also briefly discussed. The conclusions are summarized in section 4. 



2. Two- fluid cylindrical equilibria with reversed magnetic shear 

The two-fluid equilibrium states of an ideal quasineutral plasma are governed by the 
following set of equations written in Gaussian units with both 4tc and the velocity of light 
being set to unity: 

V • (nv a ) = 0, (2) 
V)v a = -VP a + (E + v Q xB), (3) 

v Q • VT a = 0, (4) 
ZiUi ^n e = n, (5) 
V x E = 0, (6) 
V • B = 0, (7) 
VxB = ^] n a q a v a = J, (8) 

where the index a denotes the particle species (a = i for ions and e for electrons); n 
is the plasma density in connection with the quasi- neutrality condition (J5|): q a is the 
charge of each particle species with Zi being the atomic number. The rest of the notation 
is standard. The energy equation (J3J), associated with the fact that for fusion plasmas 
the heat conduction along B is large and therefore the temperature becomes uniform on 
magnetic surfaces on a fast time scale, is particularly appropriate for electrons. For ions 
one alternatively can use an adiabatic energy equation: 

v 4 • VPj + 7 P*V • v, = 0. (9) 

Compared with the respective set of MHD equations (see for example Eqs. (2)- (6) of Ref. 
[TJ|) Eqs. (J2HH1) are advantageous in two respects: (i) the momentum equation includes the 
electric field and therefore the pressure gradient contribution to E can be calculated from 
this equation; this contribution is missing in the frame of MHD because E is calculated 
by Ohm's law, E + v x B = 0, and (ii) the current density J is related self-consistently 
to the fluid species velocities [Eq. (JSJ)]. 

The system under consideration is a cylindrical plasma of circular cross-section con- 
fined by a magnetic field having toroidal and poloidal components B z and Bg respectively. 
Also the velocity has toroidal and poloidal components and the electric field is radial. Be- 
cause of symmetry any equilibrium quantity depends solely on the radial distance r; 
therefore Eqs. (J2J), ©[and (0)], ©, and (JZj) are identically satisfied. Also the flow for 
both fluid species is incompressible (V • v Q = 0). Under these considerations 6 out of the 
12 scalar quantities remain free and can be prescribed. 

Adding Eq. J0J) for ions and electrons yields the MHD momentum equation 

-(p + M) + (1 _« a S = , (10) 
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where 

1/2 

I I : I I I ■ I - -!— I I I I I r 

Ma = 



UiUiivfa + n e m e v 2 ed 



B 2 

is the poloidal Mach number. Because of symmetry the toroidal velocity as well as the 
velocity shear (of both toroidal and poloidal components) do not appear in (fTU|) . It is 
convenient to use instead of © for the electrons. Therefore the slightly reduced 
set of equilibrium equations we will use in the following consists of Eqs. ©, ©, Q for 
ions only, (jSJ), ©, ©, ©, and (fTUj) . By expressing B e in terms of the safety factor, q = 
rBg/ (RqB z ), with 2irRo associated with the length of the plasma column, and introducing 
the normalized radius p = r /r® with tq corresponding to the plasma surface, Eq. (jl()j) can 
be put in the form 



P'(p) = -B,(p)B'.(p) 



1+ ( £ ^)) 2 l + [ A «") 2 + S (")- 2 ^(^) 2 W 



p ( Bz{py 

r \ q(p) 



Here, e = r /i? is the inverse aspect ratio and s(p) = (r/q)(dq/dr) the magnetic shear. 

On account of typical experimental ITB profiles we prescribe the quantities q, B z , f ie , 
Vi z and n as follows: 
Reversed magnetic shear profile: 



.(pHJi-^^ + ^M (12) 

Qc r min Qc r min 



where q c = q{r = 0), r min is the position of q m in, and Aq = q c — q m i n . The shape of the 
q profile is determined by adjusting the parameters q m i n , Ag and r min . Note that \s\ is 
proportional to Aq; therefore as Aq takes larger values the magnetic shear increases in 
both the s < and s > regions. A q profile compatible with experimental ones (see for 
example figure 10 in Ref. [Tl]) is presented in figure ^ 
Toroidal magnetic field profile: 

B z = B z0 [l + 8(l-p 2 )} 1/2 , (13) 

where B z q is the vacuum magnetic field and the parameter 5 is related to the magnetic 
properties of the plasma, i.e. for 5 < the plasma is diamagnetic. 
Gaussian-like ion poloidal velocity profile: 

v w = Av m p(l - p) exp ( - — — P ™ n ^ , (14) 

where the parameter h > is related to the velocity shear, i.e. \v' ie \ increases when h 
takes smaller values, and v^o defines the extremum of vw- 
Either peaked on axis toroidal velocity profile: 

Viz = v iz0 (l - p 3 ) 3 (15) 

or Gaussian-like Vi Z profile similar to that of (j!4)l : it is also noted that the results do not 
change if, alternative to (|TH|) . a peaked on axis toroidal velocity profile of the form 

v iz = v iz0 (l - p) exp (-J- 
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is employed; 
density profile: 

n = n (l-p 3 ) 3 . (16) 
In addition, the ion pressure can be expressed in terms of the total pressure by the relation 

Pi = XP , < A < 1 . (17) 

Since in tokamaks Mg < 0.1, the flow term in (jll|) is perturbative around the "static" 
equilibrium Mg = and therefore it can be neglected. It should be noted, however, 
that this approximation may be not good for non-circular cylindrical or axisymmetric 
plasmas because in these cases the convective term in the momentum equation depends 
on the velocity shear which in certain regions may become large (see for example the 
^-independent cylindrical and axisymmetric incompressible MHD equilibrium equations 
(23) and (22) in Refs. [l^j and [15] respectively). The following quantities then can be 
calculated self-consistently: the poloidal magnetic field, Bg = epB z /q, the magnetic shear 
s = (r/q)(dq/dr), the current density via Ampere's law, the pressure by integration of 
dUJ) and setting P(l) = 0, the ion and electron pressures Pi = \P and P e — (1 — X)P, 
the electric field by Eq. © for the ions 

E r (p) = , dPl } p) + v iz (p)Be(p) - v i6 (p)B z (p), (18) 

er n{p) dp 

its shear \E' r \ and ujexb by (JTJ. Also, the electron velocity components v ez and v e g can 
be determined by the relation J = ne(vj — v e ). It is noted here that the VPj term in 
(|18|) can be obtained in the framework of the ideal Hall-MHD model, alternatively to the 
complete two fluid one, which includes the generalized Ohm's law: 

E + vxB = — (JxB - VP e ) . (19) 

en 

Neglecting in the Hall-MHD momentum equation the convective flow term (which for the 
case under consideration corresponds to Mg = 0), the term j x B in (fT9"j) can be expressed 
in terms of the total pressure gradient: 

jxB = VP = V(P + P e ); 

then Eq. (|T9"j) leads to (fTHj). The above prescriptions and subsequent suggested calcula- 
tions consists a procedure to solve analytically the set of the two-fluid equations. The 
calculations have been performed analytically by developing a programm for symbolic 
computation ^7j in connection with [T5] . 

Inspection of ()18|) implies that in addition to the dependence of E r and E' r on the 
magnetic shear through the dPi/dp term, s is involved in the Vi Z term through the q de- 
pendence of Bg. The quantity o>exb is stronger affected by the magnetic shear because 
s is involved in both the v iz and Vig terms of uj ExB [see Eqs. (j2Ti|) and (J2T)) in section 
3]. These observations indicate that there is a synergetic contribution of magnetic shear 
and flow to E r , E' r , and c^exb- In this report results not obtainable within the framework 
of MHD will mainly be presented in next section. MHD results were reported in Ref. 
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3. Results 

We have set the following values for some of the parameters: B z o = IT, 5 = —0.0975, 
Zi = 1, r = lm, Rq = 3m, n = 5 x 10 19 part./m 3 , A = 0.6. The choice q min > 2 
was made because according to experimental evidence for q min < 2 strong MHD activity 
destroys confinement possibly due to a double tearing mode ^H]- A similar result was 
found numerically for one-dimensional cylindrical equilibria with hollow currents in [2*U] . 
Moreover in discharges with reversed magnetic shear in JET a correlation was found 
between the formation of ITBs and q m in reaching an integer value (2 or 3) [21] . The impact 
of the magnetic shear and flow profiles on the equilibrium characteristics was examined by 
varying the parameters Aq, q m i n , r min , h, vi z q and u^o m the ranges (4-14), (2-3), (0.5-0.6), 
(0.001-0.1), (10 5 -10 6 ms _1 ) and (10 4 -10 5 ms" 1 ) respectively; consequently q c = q min + Aq 
varies from 6 to 16 and it is guaranteed that M| ps M| , where = [n(mivf z +m e vl z )}/ B^, 
a scaling typical in tokamaks because B z m 10Bg and Vi Z ~ lOt^ [221 12H]- The impact 
of the variation of magnetic shear through Aq was studied by keeping r min and q min 
constant, while the impact of r min and q m i n was examined with constant Aq. 

First we will briefly report certain characteristics of the pressure and toroidal current 
density profiles which remain similar as in MHD. The total pressure profile, and therefore 
the Pi one, is peaked and for s < becomes steeper when \s\ increases as can be deduced 
from Eq. (jllj) (see also figure EJ). In addition (jllj) implies that the profile becomes steeper 
as the plasma becomes more diamagnetic, i.e. when B' z in connection with 5 in (j!3j) takes 
larger values. The J 2 -profile is hollow with its maximum located in the region where the 
q min lies as can be seen in figure El These characteristics are observed in discharges with 
ITBs an d are favorable for ITB formation. Especially for s > 2 a reversal of J z occurs 
in the s > region. This characteristic is discussed further in Ref. JT2j. It is also noted 
that a sufficient stability criterion for equilibria with reversed current density in the outer 
plasma area and monotonically increasing g-profiles was derived in Ref. |24j . 

The conclusions on the impact of the magnetic shear and flow on E r , \E' r \ and c^exb 
are reported on an individual basis in the rest of this section. 

3.1 Electric field (E r ) 

1. The electric field consists of the VP«, Vi Z and, Viq contributions in connection with 
the first, second, and third term in Eq. (|18|) . Each of these terms contributes 
of about the same order of magnitude to E r (figure EJ). This is consistent with 
experimental evidence [TT] . A similar result was obtained in a different way in [23] 
(see figure 4 therein). It is apparent from (|TS|) that E r depends linearly on f i2 and 

with the overall velocity contribution to E r , however, depending on the relative 
signs of v iz , v ie and B z . 

2. Typical E r profiles exhibit an extremum located in the neighborhood of the q m in 
position (figure EJ). 

3. Increase of \s\, by increasing Aq, makes the maximum of \E r \ to take larger values 
(figure E]). Pending on the direction (toroidal or poloidal) of the velocity and the 
shape of its profile, variation of Aq from 4 to 14 increases the values of the \E r \ 
maximum in a range that varies from 5.6% for purely poloidal flow to 48% for 
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purely peaked toroidal flow. It is reminded that, in addition to the s dependence of 
the V Pi term in (|T5|i. s contributes to E r synergetically with the Vi Z term (the Vig 
term is s independent). 

4. The larger r min the higher the values of the \E r \ maximum (for given values of Aq 
and q m in) as shown in figure |H1 Quantitatively for a variation of r min from 0.5 to 
0.6, the increase of \E r \ maximum varies from 36% to 70%. Also the position of the 
extremum (located in the vicinity of r min ) is displaced outwards. 

5. The larger q min the smaller the \E r \ maximum (figured). In particular, increase of 
q m i n from 2 to 3 (with Aq = 4 and r m j„ = 0.5), results in a decrease of the \E r \ 
maximum in the range (12%, 40%). 

6. When the flow shear increases (by decreasing h from 0.1 to 0.001) the extremum of 
E r remains practically unchanged in most of the cases considered. 

3.1 Shear of the electric field (\E' r \) 

1. As in the case of E r the contributions from the VPj-, Vi Z - and t^-related terms to 
E' r are of the same order of magnitude as shown in figure |H1 

2. The profile of E' r exhibits one local extremum on each side of the q min position 
(figure EJ). The two extrema are of opposite sign. 

3. Increase of \s\ increases both maxima of \E' r \ in most of the flow cases considered and 
this increase is larger in the region where s > than that where s < 0. This is shown 
in figure (jHJ). For certain combinations of the velocity components, however, the one 
extremum increases and the other decreases. Such a case with peaked toroidal v j 2 
and poloidal vie flow is shown in figure 

4. The larger r min the higher the maxima of \E' T \ ( figure ITTj) unless the case of poloidal 
velocity in conjunction with Gaussian-like toroidal one. Pending on the direction 
and the shape of the velocity this increase varies from from 8% to 42%. Also the 
profile of \E' r \ is displaced outwards as can be seen in figure [TTJ 

5. An increase of q m in results in a decrease of the \E' r \ extremum in the s > region in 
all of the cases considered while in the s < region this happens for V& = (figure 

H2J). 

6. By increasing the velocity shear the maxima of \E' r \ are also increased in all of the 
cases considered (figure fl~3|) . 

7. For either purely toroidal or poloidal flow, increase of the maximum absolute value 
of the velocity by a factor increases the maxima of \E' r \ by the same factor in all of 
the cases considered with the following exception: purely toroidal peaked flow for 
which the maximum in the s < region increases and the one in the s > region 
decreases. 
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8. For either purely toroidal or purely poloidal flow, inversion of the velocity direction 
causes a change in the sign of the two E' r extrema. Also, this inversion leads to (i) 
an increase of both maxima of \E' r \ for Gaussian-like Vi Z , (ii) an increase of the one 
\E' r \ maximum in the s > and a decrease of the other in the s < region and 
(iii) a decrease of both maxima for purely poloidal flow. For Gaussian-like v & the 
increase of the one |P£|-maximum in the s > region caused by inversion is greater 
than the increase of the other in the s < region. 

3.3 Shear of the E x B velocity (ujexb) 

1. The profile of uexb (Eq. QJ) possesses two maxima located the one in the s < 
and the other in the s > regions (figure EJ). Larger of the two maxima is the one 
which lies in the region of steeper pressure profile. 

2. The impact of the magnetic shear on we x b is stronger than that in MHD due to the 
VPj term of the electric field [Eq. 1)18)1] . Specifically for constant B z and arbitrary 
profiles of q, v iz and u^, Eq. (JTJ) yields at the point where E' r = 



WExB 



x (l-s)(2-s)B z pe 
enqrl(p 2 + K) 



(20) 



where 

{l-s){^-v ie ) 
WEXB - M ™ " r V + $) ( } 

The first term in ()2Uj) stems from the VPj part of E r in ([lip while the second term 
comes from the t>j 2 and Vi$ parts of E r . The subscript MHD is used to emphasize the 
similarity of ()2*T)) with the respective MHD relation derived in Ref. ^2] [equation 
(18) therein]. It is apparent the V P r related dependence of we x b on s, proportional 
to (1— s)(2— s), is stronger than the Vi Z and vw dependence proportional to 1— s; also, 
the absolute values of the VPj-, v i Z - and Vie- related terms are individually larger 
for s < than s > 0. The contribution of each of these terms to c^exb however is 
of the same order of magnitude. Note that despite of the tokamak pertinent scaling 
Vi Z ~ lOt>i0, the contributions of Vi Z and vm terms are of the same order of magnitude 
because of the factor ep/q. The "equipartition" of the three terms holds in general 
for the whole uj^xB profile obtained via the symbolic computation programme as 
shown in Fig 1151 

3. Increase of the flow via either \vi Z o\ or \vieo\ by a factor increases the maxima of 
uexb by the same factor. 

The impact of the magnetic shear through Aq and the flow on wexb is similar as that 
on E' r . Specifically: 

1. Increase of \s\ leads to larger values for the maxima of uj-^xb in most of the flow 
cases considered ( figure HE)) . There are some combinations of velocity components 
however for which the one maximum increases and the other decreases. Such a 
case is shown in figure El I n which region (s < or s < 0) the increase takes 
place depends on the particular velocity components involved and the shape of the 
toroidal velocity profile. 
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2. The larger r min the greater the WExB-maxima f figure ITSj) in the same cases as for 
E' r . The profile of cuexb is also displaced outwards as can be seen in figure fTSl 

3. Increase of q m i n causes (i) a decrease of the cjExB-maximum in the s > region in 
all of the cases considered and (ii) an increase of c^exb -maximum in the s < one 
for v ie 7^ (figure HHJ) . 

4. The larger the flow shear the greater the we x b maxima in all of the cases considered 
(figure HI- 

5- co>exb * s a ff ec t e d by the relative signs of V{g, v iz and B z as is apparent from Eqs. (J2U|) 
and (|2ip. In particular (i) inversion of the Gaussian-like toroidal velocity increases 
the maxima of cjexB (figure EH) , (ii) the maximum of co»exB m the s > region 
increases while the one in the s < region decreases due to the reversal of the 
peaked toroidal velocity and (iii) they decrease by inversion of the poloidal velocity. 
Finally for both Gaussian-like velocity components the variation caused by inversion 
is greater in the s > region. 

4. Summary and Conclusions 

In this report tokamak equilibria with reversed magnetic shear and sheared flow have 
been studied within the framework of two-fluid model in the limit of infinite aspect ratio. 
The study is based on a slightly reduced set of two-fluid equations in which the electron 
momentum equation is replaced by the respective MHD one. Neglecting the flow term in 
this equation (because in cylindrical geometry it is small for tokamaks) and prescribing 
the profiles of six free quantities in accord with ITB experimental ones, i.e. the toroidal 
magnetic field B z , the safety factor q, the toroidal and poloidal ion velocities Vi Z and Vig, 
the density n [Eqs. (fT^jl -ffTo ^) ] and the ion pressure in terms of the total pressure, Pi = 
XP, we have constructed analytic solutions in calculating self consistently the following 
quantities: P [and therefore Pi and the electron pressure P e — (1 — X)P], the current 
density and the radial electric field E r ; the electric field shear, \E' r \, and the shear of the 
E x B-velocity, c^exb [Eq. (HJ)] have also been calculated. Gaussian-like profiles for v ie 
and either Gaussian-like or peaked-on-axis ones for Vi Z have been considered. In addition, 
for reversed magnetic shear profiles the impact of s and the flow on the equilibrium 
characteristics has been examined by varying the parameters Aq which \s\ is proportional 
to, the minimum of q, q m i n , its position, r min , the extrema of the velocity components, v iz0 
and v ioo, and a parameter h which decreases with increasing velocity shear. The results 
are as follows. 

1. The pressure profiles become steeper in the region of s < 0. 

2. The profile of the toroidal current density J z is hollow and a reversal occurs in the 
outer plasma region for s > 2 in connection with appropriate values of Aq. 

3. The | E r | profile has a maximum located close to the q m in position while the \E' r \ 
and cuexb ones have two local maxima the one in the s > and the other in the 
s < regions. 
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4. The contributions associated with VPi, Vi Z , and Vi$ to E r , E' r and wexb (the VP 
contribution being missed in MHD) are of the same order of magnitude. 

5. The magnetic shear affects E r and E' r explicitly through VPi and implicitly in 
conjunction with Vi Z ; s has an additional impact on we x b m connection with Vi$. The 
explicit impact of s is stronger; in particular for B z =constant, the VPj contribution 
to ojexb at the point where E' r = is proportional to (1 — s)(2 — s) [Eq. (J2Uj) ] while 
the contribution through the flow terms is proportional to (1-s) [Eq. (121)) ]. 

6. Increase of \s\ results in an increase in the maximum of \E r \ in all of the cases 
considered. Also, the maxima of \E' r \ and o>exb increase in most of the flow cases 
considered. When either the toroidal and poloidal velocity contributions cancel each 
other or the velocity is purely toroidal peaked, the increase is greater in the s > 
region. Also pending on the direction and shape of the flow, the increase varies from 
56.4% to 323%. 

7. The larger r min the greater the maxima of \E r \, \E' r \, and cuexb- 

8. The larger q m in the smaller the maximum of \E r \ but the larger the maxima of \E' r \ 
and cuexb in the s > region. 

9. Stronger flow, by larger values of \v iz0 \ and |^o|) leads to linear increase in E r , \E' r \, 
and w ExB . 

10. The larger the flow shear (by smaller values of the parameter h) the slightly smaller 
the maximum of \E r \ but the larger the maxima of \E' r \ and uexb- 

11. E r , E' r , and co>exb are sensitive to the relative signs of v iz , v^, and B z . 

In summary, alike MHD, in the framework of two-fluid model the magnetic shear and 
sheared flow (toroidal and poloidal) act synergetically on E r , E' r , and cc'exb which may 
play a role in the formation of ITBs. However the impact of magnetic shear on these 
quantities is stronger than that in MHD due to the additional Vp contribution to the 
aforementioned terms. 
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Figure captions 

Figure ^ Safety factor profile in connection with Eq. (|T2*|) compatible with the experi- 
mental one measured in JT-60U ^3] (figure 10 therein). 

Figur^21 Pressure profiles for two values of the reversed-magnetic-shear parameter Aq 
normalized with respect to the value of P at the magnetic axis. 

Figure E3 Toroidal current density profiles for two values of Aq which show the hollow 
shape and the reversal in the outer plasma region. The profiles are normalized with 
respect to the maximum of J z for Aq = 4. 

Figure |U Profiles of the Vp, v iz and v ie contributions to the electric field, E r , showing 
that all three contributions are of the same order of magnitude. The V{ Z profile is peaked 
on axis. The profiles are normalized with respect to the extremum of VPj contribution. 

Figure EJ Increase of the normalized absolute value of the electric field extremum due to 
the variation of Aq for q min = 4, r min = 0.5, Gaussian-like v iz profile and ViQ = 0. 

Figure EH Increase of the absolute value of the E r extremum when the distance r m i n in 
connection with the q m in position becomes larger. Also, the position of the extremum is 
displaced outwards. Here, q m i n = 2, Aq = 4, and the flow is purely toroidal peaked on 
axis. The profiles are normalized with respect to the extremum for r m i n = 0.5. 

Figure [7| Decrease of the normalized \E r \ maximum when q min increases for Aq = 4, 
r m in = 0.5, peaked Vi Z and Gaussian-like localized (h = 0.001). 

Figure El Profiles of the VPj, V{ Z and contributions to the electric field shear, E' r , 
showing that all three contributions are of the same order of magnitude. The Vi Z profile 
is peaked on axis. The profiles are normalized with respect to the extremum of the VPj 
contribution in the s < region. 

Figure El Increase of the normalized | E' r | maxima caused by an increase of the magnetic 
shear in connection with variation of Aq. The plots were obtained for Vi Z peaked and 
v ie = 0. 

Figure ITU1 Increase of the \E' r \ extremum in the s > region due to the increase of Aq. 
The profiles are obtained for v% z peaked and v & 7^ and are normalized with respect to 
the value of the extremum in the s < region for Aq = 4. 

Figure HU Increase of the \E' r \ extrema as r min takes larger values for q min = 2, Aq = 4, 
and peaked purely toroidal flow. Also the positions of the extrema are displaced outwards. 
The profiles are normalized with respect to the E' r extremum for r m , n = 0.5 

Figure IT21 Decrease of the \E' r \ extremum in the s < region for purely toroidal Gaussian- 
like flow due to the increase of q m in- For this particular case of flow the variation of the 
other extremum in the s > region is negligible. 

Figure Increase of the \E' r \ extrema due to the increase of the flow shear for Gaussian- 
like Vi Z and Vie = 0. 

Figure ITU A typical c^exb profile for purely poloidal flow, normalized with respect to 
the maximum value in the s > region. 
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Figure IT31 Profiles of the VPj, Vi Z and Vi$ contributions to cuexb showing that all three 
are of the same order of magnitude. The Vi z profile is peaked on axis. The normalization 
is made with respect to the maximum value of the VPi contribution in the s < region. 

Figure ED Increase of the normalized u;e x b maxima due to the increase of the magnetic 
shear. Both velocity components have Gaussian-like profiles. 

Figure IT7I Increase of the wexb extremum in the s < and decrease in the s > regions 
due to the increase of Aq. The profiles are obtained for v iz Gaussian-like and = and 
are normalized with respect to the s < extremum of c^exb for Aq = 4. 

Figure HH1 Increase of the normalized wexb maxima due to the outward shift of the 
position of q m i n for purely toroidal Gaussian-like flow. 

Figure HHJ Increase of the c^ExB-extrema as q m i n takes larger values when both velocity 
components have Gaussian-like profiles. In this particular case the increase of the ex- 
tremum in the s > region is very small. The profiles are normalized with respect to the 
s > maximum of cc'exB for q m i n = 3. 

Figure 1201 Increase of the cuexb extrema caused by the increase of the flow shear for t> i2 
peaked and vw ^ 0. The normalization is made with respect to the c^exb extremum in 
the s > region for h = 0.1. 

Figure |^ Increase of the normalized WExB-extrema caused by inversion of a Gaussian- 
like toroidal velocity. 
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Figure 1: Safety factor profile in connection with Eq. (JT2*Jl compatible with the experi- 
mental one measured in JT-60U [Hj (figure 10 therein). 
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Figure 2: Pressure profiles for two values of the reversed- magnetic-shear parameter Aq 
normalized with respect to the value of P at the magnetic axis. 
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Figure 3: Toroidal current density profiles for two values of Aq which show the hollow 
shape and the reversal in the outer plasma region. The profiles are normalized with 
respect to the maximum of J z for Aq = 4. 




Figure 4: Profiles of the VPj, v iz and contributions to the electric field, E r , showing 
that all three contributions are of the same order of magnitude. The v iz profile is peaked 
on axis. The profiles are normalized with respect to the extremum of V Pi contribution. 



15 



E r (p)/E c 

0.425 ' 0.45 ' 0.475 ' 0.5 ' 0.525 /6 . 55 ' 0.575 \ 




y Aq = 14 



Figure 5: Increase of the normalized absolute value of the electric field extremum due to 
the variation of Aq for q min = 4, r min = 0.5, Gaussian-like v iz profile and = 0. 
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Figure 6: Increase of the absolute value of the E r extremum when the distance r min in 
connection with the q m in position becomes larger. Also, the position of the extremum is 
displaced outwards. Here, q m in — 2, Aq = 4, and the flow is purely toroidal peaked on 
axis. The profiles are normalized with respect to the extremum for r min = 0.5. 
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Figure 7: Decrease of the normalized \E r \ maximum when q min increases for Aq = 4, 
r min = 0.5, peaked v iz and Gaussian-like localized (h = 0.001). 
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Figure 8: Profiles of the VPj, Vi Z and contributions to the electric field shear, E' r , 
showing that all three contributions are of the same order of magnitude. The Vi Z profile 
is peaked on axis. The profiles are normalized with respect to the extremum of the VPj 
contribution in the s < region. 
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Figure 9: Increase of the normalized \E' r \ maxima caused by an increase of the magnetic 
shear in connection with variation of Aq. The plots were obtained for Vi Z peaked and 
v ie = 0. 




Figure 10: Increase of the \E' r \ extremum in the s > region due to the increase of Aq. 
The profiles are obtained for Vi Z peaked and ^ and are normalized with respect to 
the value of the extremum in the s < region for Aq = 4. 
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Figure 11: Increase of the \E' r \ extrema as r min takes larger values for q min = 2, Aq = 4, 
and peaked purely toroidal flow. Also the positions of the extrema are displaced outwards. 
The profiles are normalized with respect to the E' r extremum for r min = 0.5. 
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Figure 12: Decrease of the \E' r \ extremum in the s < region for purely toroidal Gaussian- 
like flow due to the increase of q m i n - For this particular case of flow the variation of the 
other extremum in the s > region is negligible. 
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Figure 13: Increase of the \E' r \ extrema due to the increase of the flow shear for Gaussian- 
like Vi Z and Vie = 0. The profiles are normalized with respect to the extremum in the 
s > region for h = 0.001. 



^exb(p)/^, 
1 



0.8 



0.6- 



0.4 - 



0.2 




0.45 



0.5 



0.55 



. 6 



Figure 14: A typical cjexb profile for purely poloidal flow, normalized with respect to the 
maximum value in the s > region. 
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Figure 15: Profiles of the VPj, Vi Z and Vi$ contributions to wexb showing that all three 
are of the same order of magnitude. The v iz profile is peaked on axis. The normalization 
is made with respect to the maximum value of the VP« contribution in the s < region. 
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Figure 16: Increase of the normalized uj-^xb maxima due to the increase of the magnetic 
shear. Both velocity components have Gaussian-like profiles. 
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Figure 17: Increase of the c^exb extremum in the s < and decrease in the s > regions 
due to the increase of Aq. The profiles are obtained for v iz Gaussian-like and v ^ = and 
are normalized with respect to the s < extremum of wexb for Aq = 4. 
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Figure 18: Increase of the normalized cjexb maxima due to the outward shift of the 
position of q m i n for purely toroidal Gaussian-like flow. 



22 



^exb(p)M 




Figure 19: Increase of the c^ExB-extrema as q min takes larger values when both velocity 
components have Gaussian-like profiles. In this particular case the increase of the ex- 
tremum in the s > region is very small. The profiles are normalized with respect to the 
s > maximum of c^exB for q m i n = 3. 
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Fi gure 20: Increase of the wexb extrema caused by the increase of the flow shear for vi z 
peaked and ^ 0. The normalization is made with respect to the c^exb extremum in 
the s > region for h = 0.1. 
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Figure 21: Increase of the normalized c^ExB-extrema caused by inversion of a Gaussian-like 
toroidal velocity. 
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